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Abstract 


la 


We explore the physics of thermal impedance matching at the interface between two dissimilar 
materials by controlling the properties of a single atomic mass or bond. The maximum thermal 
current is transmitted between the materials when we are able to decompose the entire heterostruc¬ 
ture solely in terms of primitive building blocks of the individual materials. Using this approach, 
we show that the minimum interfacial thermal resistance arises when the interfacial atomic mass 
is the arithmetic mean, while the interfacial spring constant is the harmonic mean of its neighbors. 
The contact induced broadening matrix for the local vibronic spectrum, obtained from the self¬ 
energy matrices, generalizes the concept of acoustic impedance to the nonlinear phonon dispersion 
or the short-wavelength (atomic) limit. 
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I. INTRODUCTION 


Today’s experimental techniques are opening up the possibility of tuning thermal con- 


Q. 


ductivity of materials by engineering their thermal impedance at the nanoscale 111. At these 
characteristic lengths (~10nm), thermal boundary conductance (TBC) of interfaces provide 
a major contribution to the thermal conductance of the system, making critical the under¬ 
standing of impedance matching at interfaces. Phonon transport across an interface is a 
convoluted process that involves the differing phonon modes, the short coherence lengths of 
the quantized vibrations and their broadband transport properties. Also, it involves complex 
and diverse interfacial atomic structures that depend strongly on materials and fabrication 
protocols. Several experiments {2 9| and simulations 10-15| have already shown the de¬ 
pendence of TBC with interfacial impurities, mixing, defects, chemistry or bond strength. 
Nevertheless, the standard models to calculate TBC, the acoustic mismatch model [16 1 and 

n 

the diffuse mismatch model [17[, completely neglect the properties of the interface. Al¬ 
though some work has been done to include those properties into a model 
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m a 


proper identihcation of the key physics determining TBC is still incomplete but it is crucial 
for impedance matching design at the nanoscale. This will lead the emerging held of phonon 
engineering to follow the successful steps of electronics and photonics, where engineering of 
nanoscale properties has endowed the helds with high degrees of tunability. 


While the overall goal of our study is to explore the physics of thermal impedance match¬ 
ing at interfaces covering the entire gamut from ID to 3D, from linear to non linear dispersion 
and from coherent to incoherent transport, we will start building our intuition by studying 
coherent thermal impedance matching between two dissimilar ID materials by controlling 
the properties of a single mass (Fig. [T^) or spring (Fig. [T]b) in between. This toy model 
presents a starting point to understand ballistic contributions to TBC by important factors 


already identihed in t 


bond strength 


|jy,i8 


le literature, like interfacial impurities, mixing, defects, chemistry or 


19| . In fact, some authors have used this toy model to support 


their Molecular Dynamic simulation results arguing the increase of TBC with increase of 
bond strength Q,|nl,Q. 


The results for coherent, 1-D thermal impedance matching are incredibly diverse. For 
example, to achieve maximum thermal conductance we want the interfacial impedance to 
maximize the area under the transmission function, like a broadband hlter. Following this 
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FIG. 1. a) ID interface between dissimilar materials with an arbitrary atom in between, b) 
Optimal coupling between the contacts happens when we can describe the entire heterostructure 
only in terms of building blocks of individual materials. This decomposition makes the optimal 
atomic mass the arithmetic mean of its neighbors mo = mi/2+m2/2. c) Interface with an arbitrary 
bond in between, d) Maximum thermal conductance occurs when the interfacial spring constant is 
the harmonic mean of its neighbors 1 //cq = 1 l‘lk\ +1 1 ‘lk 2 (recall that a spring half as short is twice 
as strong), which follows again from a decomposition in terms of blocks of individual materials. 


criterion, we point out that the best matching interfacial mass (mo) for the sing le mass 
junction (Fig. [T^) is the arithmetic mean between the masses of the contacts 22|. For a 
single spring junction (Fig. [Hb), Zhang et al. 20|] found that the best matching interfacial 
spring constant (/co) is the harmonic mean between the contact springs. When the goal is to 
achieve maximum phonon transmission around a hxed frequency, our expectations based on 
our knowledge of optical antireflection coatings posit that unity phonon transmission would 
require a quarter wave plate with an impedance equal to the geometric mean of its neighbors. 
In view of these diverse results, our aim here is to put these averages on a common footing 
and motivate them qualitatively in terms of the intrinsic physical properties of the junction 
itself. 


The central point of this paper is that the degree of mismatch at a single atom or bond 
interface depends on our ability to express the entire heterostructure solely in terms of prim¬ 
itive building blocks on either side. For instance, we hud that the optimal mass (Fig. [T^) 
is one that can be decomposed precisely into two half-masses arising from the materials 
on either side (Fig. [T]c). This decomposition makes the optimal mass the arithmetic mean 
of its neighbors, i.e. tuq = mi/2 -|- m2/2. For an analogous decomposition of the spring 
constant (Fig. [TJo and d), we hnd that the optimal spring constant equals the harmonic 
mean of its neighbors, i.e. l//co = l/2fci -|- l/2/c2 (Recall that springs in series add like resis¬ 
tances in parallel). Any deviation from those optimal decompositions (“Optimally Coupled 
Interfaces”-OCI) adds an extra barrier for heat carriers reducing the interfacial transmission. 
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The thermal conductance for OCI is characterized by the contact induced broadening 
matrix r(a;) extracted from the local vibronic spectrum, which generalizes the concept of 
acoustic impedance (Z). T not only includes non linear dispersion and short-wavelength 
(atomic) limit effects but its matrix character can account for the different modes or channels 
available for transport when higher dimensions are considered. Also, this character can 
include intricate chemical details at the interface, which may greatly affect the transport 

23j. It is worth emphasizing that T alone is not 


process as shown recently by Losego et at 
enough to correctly represent general phonon transport. The broadening’s Hilbert transform 
must be also included in the Green’s function to properly account for the sum rule of the 
local density of states [24 1. 


This generalization not only relates the continuum formalism 2^ with the discrete Non- 


Equilibrium Green’s Functions (NEGF) formalism 
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29l |. but also provides a way to 


extrapolate known results based on acoustic impedance to OGI. For instance, we can totally 
eliminate the interfacial reflection by impedance matching of the T matrices (more precisely, 
the projected self-energies, S), realized when T for the central layer equals the geometric 
mean of its neighbors. The generalization may also allow us to use existing techniques from 
other engineering fields in phonon engineering. For instance, broadband filter techniques 
from microwave engineering may be useful to engineer interfaces with maximum thermal 
conductance. 

The document begins by explaining the idea of splitting ID chains into primitive blocks, 
which define the properties of contacts or semi-inhnite chains (Section [IT])- Then, using the 
block concept, phonon transmission is calculated in section HTTl where it is also shown that 
maximum thermal conductivity occurs when the entire heterostructure can be expressed 
solely in terms of the building blocks on either side of the interface. These types of interfaces 
(OGI) are studied and characterized in section IIVI where it is shown that OGI generalizes 
an abrupt interface in the continuum limit with T generalizing Z. 


II. BLOCK PARTITION OF ID CHAIN AND CONTACTS 

An inhnite ID chain of masses coupled by springs (Fig. [2^) can be decomposed into 
different arrays of primitive blocks (Figs. [2]3 and [2l3 ). According to the blocks, different 
contacts, i.e. semi-inhnite chains, can be built from the same homogeneous material. As 


4 










a) ^ c) 2n\ 2n+l 2n+2 

^ ki ^ fei ! ri+l ki -^—i --^1- *■ -}- 

, ,. ; ,s 2ki SfcTj 

n ' n+l _ SLI_5zO_ 

)^WV{^)-^AAr(^)^W^r(] 


^ fci mi 
2 


|>vw<]|>wv(] 


ar=0 


il ii ll' 
2 2 ! 


FIG. 2. a) ID infinite chain of masses mi, separated by a, joined by springs with force constant ki. 
b) and c) Same chain separated into different blocks whose boundaries define different contacts, b) 
Partition into half-mass blocks with half masses across the boundary, c) Partition into half-spring 
blocks with half springs across the boundary, d) The half-spring blocks can be reinterpreted as 
half-mass blocks, provided the corresponding mass and spring constant are frequency dependent. 
The non-white background of the blocks represent this dependence. 

we will also show, we can equally view a chain of half-spring blocks as a virtual chain of 
half-mass blocks, provided that the corresponding mass and spring constant are frequency 
dependent. 

Consider a ID infinite chain of masses mi, separated by a distance a and connected by 
springs with force constant ki (Fig. [2^). Newton’s equation for the normalized displacement 
of the atom, with dimensions inverse square root of mass, i.e. [p] = 

, is given by 

OJ^miUn = -klUn-l + 2kiUn - kiUn+l- ( 1 ) 

This set of periodic equations is solved by plane waves Un = satisfying the dispersion 

relation 

cu^mi = 2ki — 2ki cos(gia). (2) 

Splitting each mass into its series equivalent mi = mi/2+mi/2, the chain can be partitioned 
into blocks with boundaries at the masses, i.e. half-mass blocks (Fig. [2]b). In this case, Eq. [ 1 ] 
and [2] can be reorganized to reflect the partition as 

2 /mi mi\ 

^ ^ + 2kiUn - kiUn+l ( 3 ) 

^ ~ cos(gia). (4) 

Note that the plane waves solving those equations represent l^piVmi propagating phonons 
of energy hw {N is the number of atoms in the chain) and carry a thermal current given by 

-phm 

J = hki sin(gia)|y4|^ = h—= fku — Vg{u:)\A\^, (5) 

^ Qj 
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where 


= 2fci sin(q'ia) (6) 

(“hm” stands for half-mass) is the non-zero entry of the broadening matrix used in NEGF 
formalism and Vg{uj) is the frequency dependent phonon group velocity. 

Similarly, we can split each spring into its series equivalent 1/ki = l/2ki + 1/2/ci, sepa¬ 
rating the chain (Fig. [2^) into blocks with boundaries at the springs, i.e. half-spring blocks 
(Fig. [2]c). The latter system is described by 


U}‘^miU2n+l = -2kiU2n + 2(^1 ki)u2n+l “ 2kiU2n+2, 


(7) 


0 — —2kiU2n-l + 2(/Ci -|- ki)u2n ~ 2/CiM2n+l- 


( 8 ) 


Solving for U 2 n and U 2 n +2 from Fq. |8]and replacing into Fq. [7]yields Fq. [T] More interestingly, 
solving for U 2 n-i and U 2 n+i from Fq. Hand replacing into Fq. |8l results in an equation similar 
to Fq. [3] 


UJ 


2 



U2n ti'U2n—2 T 2fi'U2n tlU2n+2i 


(9) 


where 


ei 


mi 



ti = 


1 


ki 



are frequency dependent coefficients and the cut off frequency is given by Ud = 2ijki/mi. 
In other words, the half-spring block chain can be interpreted as a virtual half-mass block 
chain having frequency dependent masses and springs. 

This analogy permits the extrapolation of algebraic treatments, like NFGF, from half¬ 
mass block to virtual half-mass block chains. For instance, plane waves describing the 
displacement at the boundaries of half-spring blocks U 2 n = satisfy the 

dispersion relation 

(y + y) ^ cos(gia) (10) 


and carry a thermal current 


J = htisin{qia)\A\^ = h—^\A\^ = hjjj—Vg{u:)\A\^ (11) 

2 CL 


with 

= 2ti sin(gia) 


( 12 ) 
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(“hs” stands for half-spring) the non-zero entry of the broadening matrix used in NEGF 
formalism for the virtual chain. 

Although the same inhnite chain or bulk material can be built from any block, differ¬ 
ent contacts are created from different blocks. Indeed, the block choice dehnes the edge 
of the contact, the positions in space described by displacement plane waves (block 

boundaries) and more importantly the thermal current carried by those waves. One striking 
example of the difference between half-mass and half-spring contacts arises when we connect 
them together. A phonon impinging on such an interface has non-zero probability of reflec¬ 
tion, unlike a phonon propagating in a single block chain. Note that this interface mimics a 
growth defect in a ID crystal. 

III. TRANSMISSION USING BLOCKS 

Since a set of phonons of equal energy propagating in a crystal are well represented 
by plane waves, the transmission probability of phonons impinging at an interface can be 
calculated from the ratio between the thermal currents carried by the transmitted and 
incident waves. This section presents phonon transmission calculations using the block 
concept to simplify the process. It is shown that maximum transmission at every frequency, 
and therefore maximum thermal conductance, happens when the entire heterostructure can 
be expressed solely in terms of building blocks on either side. This idea is equivalent to 
choosing the interfacial atomic mass as the arithmetic mean or the interfacial spring constant 
as the harmonic mean of its neighbors. 


A. Interface with Mass Junction 


Imagine chopping the materials of Fig. [T^ into half-mass blocks and the interfacial mass 
into a series equivalent that completes the contacts’ blocks plus some residual mass rrii, i.e. 


mi m2 

mn =- \- rrii - 

u 2 2 


(13) 


(Fig. |3h). Assuming incident, reflected and transmitted plane wave solutions, transmission 
is given by the ratio of transmitted over incident thermal currents 

Ji 
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FIG. 3. Decomposition of interfaces into blocks for transmission calculations. An upper bound 
for transmission happens when the impurity atom rrii or bond 1/ki are zero, a) Atom junction 
interface split into half-mass blocks, (b) Bond junction interface split into half-spring or (c) virtual 
half-mass blocks. In this case U 2 n = 


Un. — 


(14) 


The relationship between A and C is found by substituting the assumed solution 

Ce^q2na ^ > Q 

into Newton’s equation at the interface of Fig. [3^ (n = 0) 

2 /mi m2\ 

uj + mi + —J uo = -fciM-i -b (fci -b k2)uo - k2Ui. 

This process is simplified using Eq. 01 noting that the real part of the right hand side of 
Eq. [T5] exactly cancels a;^(mi -b m2)uo/2, which yields 


(15) 


uj'^rriiUo = i{A - B)-^ - 


hm 


(16) 


Combining this result with the fact that Uq = A + B = C, the transmission T is found to 
be a Breit-Wigner form 


T{u,mi) = 




.2 ’ 


(17) 


+ (T^™ + Tf^y 

with mi being the deviation of the interfacial mass mo from the arithmetic mean between 
the contact masses (Eq. HSj). 

Note that the largest possible transmission for every uj is obtained when m* = 0. This 
choice maximizes the thermal current flowing across the interface and the thermal conduc¬ 
tance of the system. That is, 



















In this particular case, the system becomes equivalent to an abrupt interface between con¬ 
tacts built with half-mass blocks, which is referred as an “Optimally Coupled Interface” 
(OCI). When 7 ^ 0, transmission decreases (for all cj) so rrii can be associated with an 
extra barrier lowering the thermal conductance. Specifically a delta scattering center like a 
single point impurity or defect at the interface. 


B. Interface with Spring Junction 


Imagine now chopping the contacts of Fig. [TJ3 into half-spring blocks and the interfacial 
spring into its series equivalent 


1 _ 1 1 1 

ko 2ki ki 2 / c 2 


(18) 


(Fig. ISj))- Assuming incident, reflected and transmitted plane wave solutions at the blocks’ 
boundaries, transmission is given by the ratio of transmitted over incident thermal currents 


hs 


■h F(*« 


The relationship between A and C is found by substituting the assumed solution 


j^^iqina _|_ iqina ^ Q g^^d 77. = 0” 


U2n — 


Ce^l2r 


n 


> 0 and n = 0 + 


(19) 


into Newton’s equation at the interface (n = 0 and n = O’*") for the virtual chain (Fig. [3t) 

( 20 ) 


Up'—U q- — —tiU-2 + (ti + ki)uQ- — kiUQ+ 


uj'^ —Uq+ — —kitlQ- -|- {ki t2)uQ+ — t2U2- 

This process is simplified using Eq. [101 noting that the real part of tpUQ- —U- 2 ) from Eq. 
and ^ 2 (^ 0 + — U 2 ) from Eq. | 2 T] exactly cancel uj'^eiUQ-/2 and uj‘^e 2 UQ +/2 respectively, which 
yields 

( 22 ) 


( 21 ) 


p/25 

0 = kiUQ- - kiUQ+ + i{A - B)^ 


and 


0 = —kiUQ- + kiUQj^ — iC 


-phs 

2 


(23) 
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Combining these two results with uq- = A-hB and Uo+ = C, the transmission T is given by 


r(u;, k-^) = 


4r^ 


hs'phs 
^ 2 


1 

4A? 


(r^T^") + (r^* + r^") 


2 ’ 


(24) 


with measuring the deviation of the interfacial spring ko from the harmonic mean be¬ 
tween the contact springs (Eq. [T 8 ll . 

Note that the largest possible transmission for every to is obtained when k~^ = 0. This 
choice maximizes the thermal current flowing across the interface and the thermal conduc¬ 
tance of the system. In that case, the system also becomes equivalent to an abrupt interface 
between contacts built with half-spring blocks, which is also referred as an “Optimally Cou¬ 
pled Interface”. When k^^ ^ 0, transmission (for all u) and thermal conductance decrease, 
so that k~^ can be associated with an extra barrier at the interface. 


IV. OPTIMALLY COUPLED INTERFACES 

An “Optimally Coupled Interface” (OCI) is an abrupt interface between half-mass 
or half-spring block contacts, which was proven equivalent to the single atomic or bonding 
interface with maximum possible transmission or thermal conductance (Section llllj) . This 
section shows that an OCI can be thought of as a step barrier for phonons responsible for 
the scattering due to a change in propagation medium. On the other hand, a non-OCI is 
represented by the same step barrier plus an extra barrier caused by a deviation from the 
optimal case (m* 7 ^ 0 or k~^ 7 ^ 0). The extra barrier decreases thermal conductance and 
can be associated with additional scattering mechanisms at the interface, such as impurities, 
mixing or dislocations. This section presents a useful way to visualize transmission in OCI 
from the contact broadenings and extends the concept of OCI to abrupt junctions between 
contacts built with different types of blocks. 

The section also shows that OCI generalizes an abrupt interface in the continuum limit 
without the long wavelength constraint. Moreover, the fen/A: property Z (acoustic impedance) 
is generalized by the contact property T (broadening), which unlike Z includes the atomistic 
details of the contact’s edge and the non-linear effects of phonon dispersion. This analogy 
shows a way to extrapolate previous results for interfaces in the continuum limit to the 
discrete limit by replacing abrupt interfaces with OCI and Z with T. For instance, the 
result of a thermal antireflection coating for a quarter wave length plate is obtained when 
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the plate broadening equals the geometric mean of the individual contact broadenings. 


A. Continuous vs. Discrete Limit 


The continuous medium approximation assumes that the wavelengths of interest are large 
enough (A >> a) so that the atomistic details of the media are ignorable, the dispersion is 
linear and the group velocity is constant. Within this approximation, the scattering problem 
at an interface (Table [I]-a) is solved by assuming incident, reflected and transmitted plane 
waves solutions and imposing boundary conditions on them to guarantee the validity of 
the wave equation at the interface. These conditions are nicely simplihed introducing the 
concept of acoustic impedance (Table [H-b and c [251]) 


m I H 
Z = pVg = ~ I 



= '/km. 


(25) 


From them, the ratio of the wave amplitudes is calculated and then transmission is found 
from the ratio of transmitted over incident energy currents (Table [I]-e). Scattering at these 
interfaces can be connected with medium mismatch using the reflection coefficient 


R = l-T = 


^1 — Z 2 
Z\ + Z<i 


(26) 


which vanishes only when Z\\A\/ = Z 2 \A\/ . In words, if plane waves of the same amplitude 
do not carry the same energy current in both media then some energy has to be reflected. 
That is, the scattering is solely caused by mismatch of the medium properties. 

When the long wavelength constraint is relaxed, the frequency dependent group veloc¬ 
ity, the cut-off frequency and the atomistic details of the interface affect the transmission 
(All included in calculations in section HTT]) . Nevertheless, particularizing the transmission 
calculation in section m to the optimal case {m^ = 0 for mass junction and fcj ^ = 0 for 
spring junction) displays the resemblance between an OCI and an abrupt interface in the 
continuum limit (Table [11). This suggests that OCI generalizes the continuous interface with 
T playing the role of acoustic impedance Z. Note that in the long wavelength limit T —)■ 2uZ 
(for both T^'” and T^®) and we recover the transmission result in terms of Z 


dTiTa A»a 4Z1Z2 

(0 + 0)'’ ' (Zi + Zj)^' 


(27) 
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a 


b 


c 


d 


e 

TABLE I. Parallel between an abrupt interface in the continuum limit and an Optimally Coupled 
Interface. The resembling suggests that the contact induced broadening matrix T generalizes 
the concept of acoustic impedance (Z) to the nonlinear phonon dispersion as well as the short- 
wavelength (atomic) limit. For Z, p = m/a. For T, p = m/a oi p = e/a according to the block 
choice. Note that the matrix character of T can account for all the conduction modes available in 
higher dimensions, for interactions beyond first neighbor and for tensorial properties of materials. 

Similar to the continuum limit, scattering in OCI is solely due to contact mismatch. 
That is, if plane waves of the same amplitude do not carry the same energy current in both 
contacts then some energy has to be reflected. The subtle difference from medium to contact 
reflects the fact that unlike Z, T is a contact property which ultimately depends on the block 
choice and carry information about the contact’s edge. 

After identifying contact mismatch scattering with the transmission functional dehning 
OCI (Table HJ-e), the extra term decreasing the transmission in Eq. [T7]or|2l]is associated 
with an additional source of scattering at the interface. Following this train of ideas, an 
OCI is represented by a frequency dependent step barrier for phonons responsible for con¬ 
tact mismatch scattering while a non OCI is represented by the same step barrier plus an 
extra barrier that decreases transmission and can be associated with impurities, mixing or 
dislocations (Fig. 0]). 

Note that if non-symmetric blocks are used, the simplifications in Table [l-c are not 
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Discrete 


JD —iqina 

ICJC, 



Acoustic Impedance 
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FIG. 4. a) Optimally Coupled Interface and its representation as a single barrier for phonons, which 
is responsible for the scattering due to propagating waves changing medium, b) Non OCI and its 
representation as a step barrier plus an extra barrier caused by a deviation from the optimal case. 
This extra barrier decreases thermal conductance and can be associated with impurities, mixing 
or defects at the interface. 


possible for every u. That is, abrupt interfaces between contacts built with non-symmetric 
blocks do not resemble the equations in the continuum limit. 


B. Characteristics of OCI 


OCI transmission can be visualized from the contact broadenings using the reflection 
coefficient, 


R = l-T 


hi-Ta 
hi + Ts 


(28) 


Unity transmission (T(a;*) = 1) is obtained when T’s match (ri(a;*) = r 2 (a;*)) at a particular 
frequency ca*. That is, a phonon with energy hcu* does not see the interface. Null transmission 
{T{u) = 0) is obtained if any of the T’s becomes 0 or imaginary. This defines a cnt off 
freqnency Uc over which phonons do not propagate in the contact. Also, note that the 
freqnency dependence of T disallows the possibility of matching different contacts at every 
oj, making scattering nnavoidable. 

When a contact is built with half-mass blocks, Ti from Eq. [6] can be rewritten using the 
dispersion relation as 


T 


hm 

1 




(29) 


with cjci = 2y^ki/mi. This concave function vanishes at cn = 0 and oo = cja and has a 
maximum value 2ki at frequency (Fig.Eti). With this function in mind, transmission 

of an OCI between half-mass block contacts (Fig. [S]!) can be visnalized from a plot of the 
real part of both contact broadenings (Fig. [5^) and Eq. [281 Figs. [5^ and|5]i show the case 
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FIG. 5. Transmission of OCI visualized from the broadening of the contacts (Ti and r 2 ). a), b) and 
c) F’s for half-mass contacts, half-spring contacts and half-spring, half-mass contacts respectively. 
The dotted lines represent r 2 . d) Characteristic transmission function for all the cases before. At 
the particular frequency cj*, where Ts intersect each other, transmission becomes unity and the 
materials match. 


when the contacts match at a particular frequency ca*, i.e. when T{u^) = 1. This is only 
possible if T’s intersect, which requires ki < ^2 and ojd > ujc 2 or ki > k 2 and Ud < ujc 2 - The 
intersection frequency cn* is found by equating ri(ci;*) = r 2 (ci;*) as 


4(m.ifci—m.2A:2) 


if mi 7 ^ m2 


OJ 


2 

* 


< [ 0 , min(a;ci, <^c 2 )] T mi = m 2 and ki = /c 2 
never if mi = m 2 and ki 7 ^ /c 2 


Note that when acoustic impedances match the intersection frequency is zero, which makes 
sense since in the continuous regime any wave looks like a zero frequency wave. When the 
contacts do not match or T’s never intersect, the transmission never reaches unity. This 
happens if ki > /c 2 and cUd > or ki < k 2 and cOd < uJc 2 - Another interesting case occurs 
when the cut-off frequencies are equal, which makes Ti oc r 2 and therefore transmission is 
constant. 

When a contact is built with half-springs blocks, T 1 from Eq. [T2] can be rewritten using 
the dispertion relation as 


= Akd 




/ i2 / ,2 

- U 


(30) 


This convex function vanishes at cn = 0, is cx) at a; = and has a slope at a; = 0 of 2^/klm[ 
(Fig. [5 ]d). Similar to the former case. Fig. |5 ]d and d show the case when the contacts match. 
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which is only possible if kirrii > ^ 2^12 and ojd > 0^02 or kimi < k 2 m 2 and ojd < 0 Jc 2 - The 
intersection frequency cu* is found equating ri(a;*) = r 2 (a;*) as 


u 


2 

* 


^ kik2 (miki-m2k2) 

17111712 k'f — k2 

< [ 0 ,min(a;ci,a;c 2 )] 

never 

v. 


if ki ^ k2 

if ki = k 2 and mi = m 2 
if ki = ^2 and mi 7 ^ m 2 


The contacts do not match if kimi < /c 2 m 2 and Ud > ujc 2 or kimi > /c 2 m 2 and ood < 0 Oc 2 
and transmission is constant if Ud = i^c 2 - 


C. Other OCI 

Since the key to obtain an OCI is the use of symmetric blocks, one can imagine that an 
abrupt interface between contacts built with different types of block is also an OCI. In fact 
this is shown starting from a more general interface where two parameters mo and ko can be 
varied (Fig. E]). Using the block concept to define the contacts and the impurities, phonon 
transmission is found to be 

ko mo 


-4 -2 0“ II 0+ 1 2 



2 2 

FIG. 6. ID interface between dissimilar materials with an arbitrary bond and atom in between. 
An OCI between half-spring and half-mass contacts arises when m* = 0 and 1/ki = 0. However, 
this OCI is an upper bound only for thermal conductance of systems in which mo vary arbitrarily 
and 1/ki = 0 or reciprocally. When ko and mo can vary together, interferences enter to the picture 
and the upper bound is lost. Interface with varying mass mo = m* ^ and spring ^ = -^ + 

The ID chain modeled is obtained by combining the springs and masses when possible. 


T 


YhsY 


hm 

2 


9 , 1 rf® ^ 


+ 




(31) 


From this system, two single junction interfaces can be defined by setting k~^ (or m*) to 
zero and letting mo (or ko) vary. The cancelations in the denominator of Eq. [31] expose an 
OCI when m* = 0 (or k~^ = 0) which is an upper bound of thermal conductance for any 
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FIG. 7. System consisting of two mediums sandwiching a third one characterized by the broad- 
enings Fi, r 2 and Fq respectively. Similar to the antireflection coating condition, transmission is 
unity when L = A/4 and Fq = -v/FiT^- 

other choice of m, (or k^^). Similar to the cases in the last subsection, plotting Fj'® and 
reflects some transmission characteristics (Fig. [St; and[5li). If the system is studied as a 
whole and both impurities do not vanish, the transmission of the OCI is not necessarily an 
upper bound for every oj. This does not contradict the previous definition of OCI because 
the interface allows two parameters to vary. 
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D. Extrapolation Example: Beyond Single Interfaces 


In subsection IIV Al we showed that OCIs generalize continuous interfaces to the discrete 
limit, which allows the extrapolation of known results between limits by changing interfaces 
with OCI and acoustic impedance Z with broadening F. This analogy may provide a way 
to endow phonon engineering with existing design criteria from other engineering helds. For 
instance, broadband hlter techniques from microwave engineering may be useful to engineer 
interfaces with maximum thermal conductance. As an example of the generalization, let us 
consider a system consisting of two mediums sandwiching a third one with impedances Zi, 
Z 2 and Zf) respectively. Recall that reflection is eliminated when the coupling medium has 
length of a quarter wavelen gth and impedance Zq = yJZxZ^- Applying the extrapolation 
rules to the known solution 2^, the system turns into Fig. [7] and its transmission is given 


by (the result can also be obtained following the process in section 


T = 


4 G 


cos2(goT) + ^j sin2(goT) 


a + l 


(32) 


with F’s defined according to the block choice in each particular region. Similar to the 
antireflection condition, transmission is unity when L = A/4 and Fq = vTiF^. 

Unlike the impedance formalism, i.e. Eq [32] replacing F with Z, the broadening formalism 
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FIG. 8. Comparison of the transmission function predicted by Eq. [32] (dotted line) and its 
counterpart with the long wavelength constraint, i.e. Eq [32] replacing F with Z (solid line). 
The system consists only of half-mass blocks (Fig. [7|) with nii = 10“^® kg, m2 = 10“^^ kg, 
ki = k2 = 10 N/m and N}, = 3. mo and k^ were chosen to guarantee T = l@a; = 4x 10^^ rad/s, 
i.e. mo = 1.1779X 10“^® kg and k^ = 7.0338 N/m. The line with crosses represents the transmission 
of a system consisting only of half-spring blocks with the same parameters. 

(Eq. I3^ includes the effects of non linear dispersion and atomistic details. A comparison of 
the transmission predicted by both formalisms is shown in Fig. (8] Note that as the frequency 
increases and the non linearity of the dispersion becomes important, the transmission func¬ 
tions separate from each other. Also note that different atomistic details at the interface, 
dehned by our block choice, generate different transmission functions. 


V. CONCLUSION 

This paper showed that the degree of mismatch at a single atom or bond interface depends 
on our ability to express the entire system solely in terms of building blocks on either side. 
Based on this concept, we argued that maximum thermal conductance happens when the 
mass or bond junctions are the arithmetic or harmonic mean of its neighbors respectively. 
Any deviation from those Optimally Coupled Interfaces (OCI) adds an extra barrier for heat 
carriers reducing the interfacial transmission and thermal conductance. We also showed that 
OCI and contact broadening (F) generalize continuous interfaces and acoustic impedance 
{Z) to the nonlinear phonon dispersion as well as the short-wavelength (atomic) limit. This 
generalization not only relates the continuum formalism with the discrete NEGF formalism, 
but also provides a way to extrapolate previous results based on acoustic impedance to 
OCI. This may allow us to use existing techniques from other engineering helds to phonon 
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engineering. For instance, broadband filter techniques from microwave engineering may be 
useful to engineer interfaces with maximum thermal conductance. 
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